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ABSTRACT 

Very deep near-infrared observations in the Faint InfraRed Extragalactic Survey (FIRES) have recently 
uncovered a significant population of red galaxies at redshifts z > 2. These distant red galaxies (DRGs) are 
efficiently selected by the criterion 7, — /T, > 2.3. We use Chandra data to examine the X-ray emission from DRGs 
in the 5' x 5' FIRES MS 1054-03 field. Two of 40 DRGs with K, < 22 ai-e detected by Chandra, and we infer that 
5^2 % of DRGs host active nuclei with Lx > l-2x 10"*^ ergs s^' . This fraction is smaller than that inferred from 
optical and near-IR spectroscopy, probably largely due to strong spectroscopic selection biases. By stacking all 
undetected DRGs we find that their average X-ray flux in the 0.5 — 8keV band is w 4.6 x 10^'^ ergs s^' cm~^. 
The detection is only significant in the soft (0.5 — 2keV; 3.4a) and full (0.5 — 8keV; 3.2a) energy bands. The 
mean detection may result from star formation, the presence of low luminosity AGN, or a combination of both. 
Assuming the detection is due exclusively to star formation, we find an average instantaneous star formation 
rate of 214 ± 68 (random) ±73 (systematic) Mq yr^', in excellent agreement with previous results from spectral 
energy distribution fitting when constant star formation histories are assumed. These results may imply that DRGs 
contribute significantly to the cosmic star formation rate at z « 2.5. However, the mean X-ray flux strictly provides 
only an upper limit to the star formation rate due to the uncertain contribution of low luminosity, possibly obscured 
AGNs. Observations at other wavelengths are needed to provide independent estimates of the star formation rate 
of DRGs. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

Recent studies have demonstrated that galaxies at z > 2 can 
be efficiently selected by the simple observed near-infrared 
] (NIR) color criterion Js — Ks > 2.3 (Franx et al. 2003; van 

■ Dokkum et al. 2003). These distant red galaxies (DRGs) 
are typically very faint in the rest-frame ultra-violet (UV), and 
are complementary to the UV-selected Lyman break galaxies 
(LBGs). Studies of their broad-band spectral energy distribu- 

[ tions (SEDs) (Forster Schreiber et al. 2004a) and their rest- 
. frame optical emission lines (van Dokkum et al. 2004) indicate 
' that at a given rest-frame optical luminosity DRGs are dustier, 
more massive, and have higher ages than LBGs. Their star 

■ formation rates are quite uncertain, as they depend on the as- 
' sumed star formation histories: median values for DRGs with 

A" < 21.7 range from ~ 170Mq yr~^ for constant star forma- 
tion histories to ^ 23 Mq yr^ ' for declining models (Forster 
Schreiber et al. 2004a). 

In this paper we study the X-ray emission from DRGs, using 
Chandra observations of the field centered on the foreground 
cluster MS 1054-03 at z = 0.83. This field is one of the two 
Faint InfraRed Extragalactic Survey (FIRES) fields and its pop- 
ulation of red z > 2 galaxies has been described by van Dokkum 
et al. (2003; 2004) and Forster Schreiber et al. (2004a). The 
X-ray properties provide information on the fraction of Ac- 
tive Galactic Nuclei (AGN) among DRGs. Spectroscopy in- 
dicates that this fraction could be substantial (van Dokkum et 
al. 2003, 2004), but this interpretation is uncertain because 
current spectroscopic samples are very incomplete. Further- 
more, the stacked X-ray flux of undetected sources provides 



an additional constraint on the average instantaneous star for- 
mation rate. Where needed we assume fi,,, — 0.3, f^A = 0.7, 
and Ho = 70kms^' Mpc^'. Galaxy identifications refer to the 
Forster Schreiber et al. (2004b) FIRES catalog of the MS 1054- 
03 field, which includes photometric redshifts for all sources. 

2. OBSERVATIONS AND SAMPLE SELECTION 

The X-ray data were obtained from an archival 9 1 ks expo- 
sure with the Chandra ACIS-S3 detector. The data were ob- 
tained by the ACIS GTO team to study the foreground cluster 
(Jeltema et al. 2001). The reduction of the data and an analysis 
of the point sources in the field are presented in Johnson, Best, 
& Almaini (2003). After removal of periods with strong back- 
ground flaring the effective exposure time is 74 ks. Images of 
the full energy band (0.5-8 keV) were made, as well as of the 
soft (0.5-2 keV) and hard (2-8 keV) bands. 

The DRGs were selected from the 5' x 5' area imaged with 
the VLT Infrared Spectrograph And Array Camera in the Js, H 
and Ks bands as part of the FIRES project. HST WFPC2 data 
is available in the Veoe and /814 bands, as well as VLT FORSl 
data in the Bessel U, B and V bands. The reduction, source 
detection, and photometry are described in Forster Schreiber et 
al. (2004b). There are 45 sources with 7, - > 2.3, < 22, 
and photometric weights (i.e., the fraction of the maximal expo- 
sure time spent on-source) in J, and Ks greater than 0.1. Figure 
1 shows the selection limits. From this group, we removed 2 
sources because they have photometric redshifts z < 1.8, and 
three sources within ^ 1' of the center of the cluster, to mini- 
mize effects of the diffuse cluster light. The final sample thus 
consists of 40 DRGs. 
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Fig. 1. — Color-magnitude diagram for all objects in the VLT field. 
Chandra point sources from Johnson et al. (2003) are shown as sohd 
circles. The open circle shows object 313, which is a A.la detection. 
Our selection criteria are shown with lines. The fraction of Chandra 
sources among galaxies with Js—Ks> 2.3 and < 22 is 5 %. 



3. DIRECT DETECTIONS: AGN FRACTION 

Johnson et al. (2003) selected point-sources in the MS 1054- 
03 field using WAVDETECT with a rather conservative thresh- 
old (giving ~ 1 spurious source over the entire field), and we 
first determined the overlap between DRGs at Ks < 22 and X- 
ray sources in the Johnson et al. catalog (large solid circles in 
Fig. 1). Only one DRG (object 1100; Zphot = 2.9; K, = 20.7) 
was found to contain an X-ray source identified by Johnson et 
al. (object 30 in their catalog). 

Next we examined the prevalence of X-ray sources fainter 
than the conservative limits imposed by Johnson et al. The 
X-ray fluxes were determined using apertures of 4" diameter, 
which contain ^ 70 % of the flux for point-sources over our 

entire 5' x 5' FIRES field ^. For each DRG a local background 
was determined from the mean flux in 30 randomly placed aper- 
tures near the location of the DRG. The locations of these aper- 
tures were constrained to be non-overlapping, at a distance of 
7" — 20" from the DRG, and at least 4.5" from point sources 
identified by Johnson et al. In addition to galaxy 1100 we find 
one other object (313; Zphot = 2.0; = 20.2) with a significant 
(w 4.7ct) X-ray flux. 

The SEDs of both objects are shown in Fig. 2, along with 
the median SED of all other DRGs (in the observed frame; for 
rest- frame SEDs see Forster Schreiber et al. 2004a). The SED 
of object 313 is very similar to the median, suggesting that the 
AGN does not contribute significantly to the continuum flux. 
However, galaxy 11 00 has much bluer optical-NlR colors than 
the majority of DRGs indicating that its active nucleus may ef- 
fect the SED, in particular in the rest-frame ultra-violet. 

The 4fT limit of 5.3 counts (after background subtraction) 
corresponds to an X-ray flux of 3.7 x 10~'* ergs s ~ ^ cm~^ (see 
below), and a K-corrected luminosity of 1.2 x 10"*^ ergs s~^ in 

^ http://cxc.harvard.edu/cal/Acis/Cal_prods/psf/Memo/s8.html 



the rest-frame 2 — lOkeV band for the median (zphot) = 2.4 of 
the DRGs in our sample. We conclude that the fraction of AGN 
withLx > 1.2x lO'^^ergss"^ among DRGs is 5^2% (assuming 
Poisson statistics). 

We note here that both spectroscopically-identified AGNs in 
the MS 1054-03 field have dropped out of the DRG sample af- 
ter recalibration of their photometry (see van Dokkum et al. 
2004). Interestingly only one of the two (object 1035, with 
Js —Ks = 2.27; labeled in Fig. 1) is detected with Chandra, sug- 
gesting that a multi-wavelength approach is necessary to accu- 
rately determine the AGN fraction among red galaxies at z > 2. 
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Fig. 2. — Spectral energy distributions of the two DRGs associ- 
ated with X-ray sources of > 1.2 x 10*^ ergss^'. The solid line 
shows the median SED of all DRGs; dotted lines show the 25- and 75- 
percentiles of the distribution. The SED of 313 (open points) is very 
similar to those of other DRGs, but the SED of 1 100 (filled points) is 
likely affected by emission from the nucleus, in particular in the UV. 

4. STACKED EMISSION OF NON-DETECTIONS 

The low X-ray background allows very efficient stacking 
analyses of sources that individually are too faint to be detected 
(e.g., Alexander et al. 2003; Brandt et al. 2001; Malhotra et al. 
2003; Reddy & Steidel 2004). We averaged the 38 flux mea- 
surements of the undetected DRGs, effectively increasing the 
exposure time to 2.8 Ms. We also averaged each of the 30 sets 
of background measurements. These 30 random stacks have 
the same properties (background and point spread function) as 
the red galaxy stack, and are used to assess the significance of 
the flux in the averaged red galaxy apertures. 

The results are listed in Table 1. The average background 
and significance in each band was determined in the following 
way. Each of the 30 sets of random apertures was averaged in 
the same way as the DRG apertures, providing 30 independent 
measurements of the mean of 38 random apertures. The mean 
of these 30 measurements gives the background value, and their 
rms is used to calculate the significance of the mean detection. 
After background subtraction, the mean emission from DRGs 
is 0.54 counts (3.2a- significance) in the full band, 0.41 counts 
(3.4cr) in the soft band, and 0.14 counts (Ict) in the hard band. 
Images of the average detections are shown in Fig. 3. 

A concern is that the mean flux is dominated by one or two 
luminous sources which are just below the detection threshold. 
Figure 4 shows the distribution of the X-ray fluxes of our un- 
detected DRG sample after background subtraction, as well as 
the average distribution of the 30 random samples. The DRGs 
show a tail of positive fluctuations and a deficit of negative fluc- 
tuations, consistent with the hypothesis that most sources con- 
tribute to the mean detection. 
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Fig. 3.— Stacked images of DRGs in the full (0.5 - 8keV), soft 
(0.5 — 2 keV) and hard (2 — 8 keV) bands. The size of each image is 
10" X 10". The DRGs are detected in the soft (3.4 a) and full bands 
(3.2 a). The detection in the hard band is not significant. 




counts 

Fig. 4. — Distribution of X-ray emission among red galaxies (solid 
line) and the average histogram for the random samples (broken line), 
after background subtraction. Compared to the random apertures the 
DRGs show a deficit of negative fluctuations and an excess of positive 
fluctuations. 

We empirically test how much flux is lost by our adopted 
4" apertures by repeating our measurements of the mean X- 
ray emission using differently sized apertures. The results are 
listed in Table 1 . As expected, the flux is a function of aperture 
size. Based on the values in Table 1 and the ACIS documen- 
tation we estimate that the 4" apertures miss 30 ± 10 % of the 
flux at the average location of the DRGs (~ 2'.5 from the op- 
tical axis). This aperture effect is partially countered by weak 
lensing of the background DRGs by the foreground z = 0.83 
cluster. This effect is « 15 % on average (see Forster Schreiber 
et al. 2004a). The aperture- and lensing-corrected mean flux of 
the DRGs is 0.66 ±0.21 (random) ±0.11 (systematic) counts 
in the full band. The systematic error reflects the uncertainties 
in the two corrections. 

The average X-ray detection of the DRGs provides a con- 
straint on their instantaneous star formation rate. In general. 
X-ray emission of galaxies is due to a combination of hot gas 
(e.g., from supernova bubbles and winds), low mass X-ray bi- 
naries, high mass X-ray binaries (HMXBs), and possibly an ac- 
tive nucleus. HMXBs are short-lived as they trace the prompt 
formation rate of massive stars with ~ SM©, and it ap- 
pears that for galaxies with star formation rates in excess of 
~ 5Mq yr~^ they dominate the X-ray luminosity, particularly 



above 2keV (Grimm et al. 2003). Several recent studies have 
therefore explored the relation between the instantaneous star 
formation rate and the rest- frame 2-lOkeV X-ray luminosity 
(e.g., Ranalli et al. 2003, Grimm et al. 2003). Here we use flie 
Grimm et al. (2003) relation, SFR« 1.49 x lO""*" x Li-iokev- 
The Grimm et al. calibration gives ^ 30 % lower star formation 
rates than the Ranalli et al. calibration. 

We convert the observed counts in the 0.5 — 8keV energy 
band to a rest-frame 2 — lOkeV luminosity in the following 
way. First we determined the average conversion from counts 
to flux from the sources hsted in Johnson et al. (2003). The con- 
version varies by less than 10 % over the ACIS field. Their con- 
version assumes a spectrum with photon index F = 1 .7, and we 
use the Portable Interactive Multi-Mission Simulator (PIMMS) 
to calculate the effect of changing F from 1.7 to 2.0, which 
is the appropriate value for the > 2keV spectrum of HMXBs. 
Applying this Ri 15 % correction to the empiricaUy deter- 
mined conversion factor we obtain Fx(0.5 — SkeV) = 6.97 x 
counts (0.5 — SkeV), with Fx in units of 10~'^ergss~' cm~^. 
The lensing- and aperture-corrected mean ORG flux is then 
4.6 x 10"'^ ergs s"' cm"^ in die 0.5-8 keV band. Next, the flux 
is K-corrected to the rest-frame 2-10 keV band (again assuming 
F = 2.0), using z = 2.4. Converting the rest-frame flux to a rest- 
frame luminosity (for z = 2 .4) we find = 1 -4 x 1 0^^ ergs s ^ ' , 
and the Grimm et al. (2003) calibration gives an instantaneous 
star formation rate of 214 ± 68 (random) ±73 (systematic) 
M0yr~'. The systematic error is due to uncertainties in the 
aperture- and lensing correction, uncertainties in the conversion 
of X-ray luminosity to star formation rate, and the scatter in the 
relation between X-ray luminosity and star formation rate. It 
was assumed that these uncertainties can be added in quadra- 
ture. We note that PIMMS predicts a ~ 4 : 1 ratio for the soft 
vs. hard band counts, consistent with our non-detection in the 
hard band. 



5. DISCUSSION 

Our estimate of the star formation rate in DRGs based on X- 
ray emission can be compared to previous estimates based on 
SED fitting and emission lines (van Dokkum et al. 2004; Forster 
Schreiber et al. 2004a). The main uncertainty in the SED 
fits is the assumed star formation history. Forster Schreiber 
et al. (2004a) find a median (mean) instantaneous star for- 
mation rate of 170 (190)MQyr-i for DRGs with K, < 21.7 
in the MS 1054-03 field assuming continuous star formation 
models, and only 23 (69)MQyr~' for declining models with 
T = 300 Myr. As shown in van Dokkum et al. (2004) the in- 
clusion of Ha measurements does not break the degeneracy be- 
tween star formation history, dust, and instantaneous star for- 
mation rate. Our result is in excellent agreement with the high 
star formation rates derived from continuous (dusty) models, 
and may imply that DRGs contribute significantly to the cos- 
mic star formation rate at z ss 2.5. 

It is also interesting to compare the mean X-ray lumi- 
nosity of DRGs with that of Lyman break galaxies (LBGs). 
Reddy & Steidel (2004) find a mean X-ray flux of - 5 x 
10~"^ergss~' cm~^ in the 0.5 — 2keV band and an imphed 
SFR of 47 Mq yr-i for LBGs at 2 < z < 2.5 in the Chandra 
Deep Field North. As the redshift range of these LBGs and our 
DRGs is roughly similar we can compare the fluxes directly, 
limiting systematic uncertainties. The mean flux of the DRGs 
is 2.3 X 10~'^ ergs s~' cm~^ in the 0.5 — 2keV band, a factor 
^ 5 higher than that of LBGs. The interpretation is unclear, as 
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the two samples have not been matched in observed K magni- 
tude, their luminosity functions may be different, and the space 
density of DRGs is still uncertain. 

Our estimate does not include one significant systematic un- 
certainty, likely to be present in most attempts to infer high red- 
shift star formation rates from X-ray luminosities. For nearby 
galaxies the angular resolution of Chandra can be used to sub- 
tract out the possible contamination from a central active nu- 
cleus, but this is not possible for high redshift objects. All the 
X-ray luminosity attributed to HMXBs could in fact come from 
a population of low luminosity AGN. Since we do not have 
optical spectra for most of our obects and we do not have the 
X-ray counts and count rates to impose constraints based on X- 
ray spectral and variability information, we cannot directly rule 
this possibility out. 

As a general caution, we remark that even if we did have opti- 
cal spectra, it might be quite difficult to rule out the presence of 
a low luminosity AGN. If we assume that the relation between 
the mass of a galaxy bulge and the mass of the bulge's cen- 
tral black hole (Ferrarese & Merritt 2000; Gebhardt et al. 2000) 
is already estabhshed for these objects, then the red galaxies in 
our sample have typical black hole masses ~ 10^ M0 . The mean 
X-ray luminosity of ^ 10"*^ ergs s^' is a factor ^100 below 
the typical X-ray luminosity for unobscured AGN with similar 
black hole masses, i.e., the nuclei of these galaxies are heavily 
obscured or the possible black holes in these galaxes are ac- 
creting at low, sub-Eddington rates. In either case, if one does 
not have the angular resolution to separate out the nucleus, the 
starlight from the bulge can completely dominate the AGN hght 
(e.g., Moran et al. 2002). The lack of obvious AGN features in 



a spectrum is thus not conclusive evidence against an AGN con- 
tribution if the aperture used contains most of the galaxy and 
the signal to noise ratio of the spectrum is low. A better (but 
still circumstantial) argument that the X-ray luminosity repre- 
sents star formation might instead be the concordance of several 
independent star formation indicators, e.g., between the radio, 
de-reddened UV, and X-ray in Reddy & Steidel (2004). Further 
independent constraints on the star formation rate in DRGs are 
thus urgently needed. 

Studies of larger samples will provide a better measurements 
of the number of luminous AGNs: whereas our initial spec- 
troscopy indicates AGN fractions of 30 — 50 % (van Dokkum et 
al. 2003, 2004) the Chandra data suggest that this is largely due 
to the strong spectroscopic selection bias toward galaxies with 
bright emission lines. It is also important to directly assess Ihe 
effect of AGNs on the broad-band SEDs (and hence estimates 
of star formation rates, ages, and dust content). Planned obser- 
vations with NICMOS and ACS on HST of MS 1054-03 will 
provide constraints on the point-source contribution to the con- 
tinuum from the rest-frame UV to the optical. Finally, deeper 
X-ray observations are needed to further constrain the hardness 
ratio of the mean detection and to detect more individual X-ray 
sources; our results and those of Barger et al. (2003) suggest 
that in 1-2 Ms exposures many of the z > 2 red galaxies may be 
(faint) Chandra sources. 
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TABLE 1 
Stacked X-ray Emission of DRGs 



Band 




Source 


2" 


3" 


4" 


6" 


Soft 




J-K 


0.36 


0.76 


1.13 


2.05 


(0.5- 


■2keV) 


Backgr. 


0.18 


0.41 


0.72 


1.62 






RMS 


0.07 


0.11 


0.12 


0.19 


Hard 




J-K 


0.33 


0.72 


1.18 


2.55 


(2-8 


keV) 


Backgr. 


0.26 


0.61 


1.04 


2.37 






RMS 


0.06 


0.09 


0.14 


0.22 


Full 




J-K 


0.68 


1.48 


2.30 


4.60 


(0.5- 


-8keV) 


Backgr. 


0.43 


1.02 


1.76 


3.99 






RMS 


0.09 


0.14 


0.17 


0.28 



